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Abstract 

In this paper we construct a self-consistent renormalization group procedure 
for MHD turbulence in which small wavenumber modes are averaged out, 
and effective mean magnetic field at large wavenumbers is obtained. In this 
scheme the mean magnetic field scales as /c -1 / 3 , while the energy spectrum 
scales as /c~ 5//3 similar to that in fluid turbulence. We also deduce from the 
formalism that the magnitude of cascade rate decreases as the strength of the 
mean magnetic field is increased. 
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I. INTRODUCTION 



Kolmogorov hypothesized that the energy spectrum E(k) of fluid turbulence in the in- 
ertial range is isotropic and is a power law with a spectral index of —5/3, i.e., 

E(k) = K Ko U 2 ^k~ 5 / 3 (1) 

where Kk is an universal constant called Kolmogorov's constant, k is the wavenumber, 
and IT is the nonlinear energy cascade rate. Note that II is equal to the dissipation rate 
and also the energy supply rate of the fluid. Experiments , simulations 0] , and some of 
the analytical calculations based on Direct interaction approximation fUJ, renormalization 
group (RG) techniques |p HTl|1 , self-consistent mode coupling etc. are in good agreement 
with the above phenomenology. 

In this paper we will discuss the energy spectrum in magnetohydrodynamic (MHD) 
turbulence. In MHD there are two fields, the velocity field u and the magnetic field B = Bo + 
b, where B is the mean magnetic field or the magnetic field of the large eddies, and b is 
the magnetic field fluctuation. One usually uses Elsiisser variables z ± = u ± b. Here the 
magnetic field has been written in velocity units (6/ y/^jfp, where p is the density of the 
fluid). We also assume that the plasma is incompressible. 

There are two time-scales in magnetofluid: (i) nonlinear time-scale l/(kz^) (similar to 
that in fluid turbulence) and (ii) Alfven time-scale 1/ (kBo). Kraichnan |13[ and Dobrowolny 



et al. |Tj| argued that the interacting z£ and z^ modes will get separated in one Alfven 
time-scale because of the mean magnetic field. Therefore, they chose Alfven time scale 
t a — {kBo)" 1 as the relevant time-scale and found that 

n+ « ir « ^-E + (k)E-(k)k z = n. (2) 

where are the cascade rates of z£. If E + {k) ~ E~(k), then the above equation implies 
that 

E + (k) ss E~(k) S3 (5 n) 1/2 &T 3/2 (3) 



In absence of mean magnetic field, the magnetic field of the largest eddy was taken as 
Bq. Kraichnan |TJ| also argued that the fluid and magnetic energies are equipartitioned. 
The above phenomenology is referred to as Dobrowolny et al.'s generalized Kraichnan (KD) 
phenomenology. 

If the nonlinear time-scale r^ L rs kz^ is chosen as the interaction time-scales for the 
eddies z£, we obtain 

n^^) 2 ^)^, (4) 

which in turn leads to 

E ± {k) = ir ± (n ± ) 4 / 3 (n^)- 2 / 3 A;- 5 / 3 , (5) 

where are constants, which we will refer to as Kolmogorov's constants for MHD tur- 
bulence. Because of its similarity with Kolmogorov's fluid turbulence phenomenology, this 
phenomenology is referred to as Kolmogorov-like MHD turbulence phenomenology. This 
phenomenology was first given by Marsch [|Tj|, Matthaeus and Zhou ||16|| , and Zhou and 



Matthaeus |T7| (it is a limiting case of a more generalized phenomenology constructed by 
Matthaeus and Zhou |Tj|, and Zhou and Matthaeus [|IjJ]). It is implicit in these phenomeno- 



logical arguments that KD phenomenology is expected to hold when Bq ^> JkE ± (k), while 



Kolmogorov-like phenomenology is expected to be applicable when Bq <C ^kE ± {k). 

In the solar wind, which is a good testing ground for MHD turbulence theories, Matthaeus 



and Goldstein [18] found that the exponent of the total energy is 1.69 ± 0.08, whereas the 



exponent of the magnetic energy is 1.73 ± 0.08, somewhat closer to 5/3 than 3/2. This is 



more surprising because Bq ^> ykE ± (k) for inertial range wavenumbers in the solar wind. 
The numerical simulations also tend to indicate that the Kolmogorov-like phenomenology, 
rather than KD phenomenology, is probably applicable in MHD turbulence |fT9"] . Hence, the 
comparison of the solar wind observations and simulation results with the phenomenological 
predictions appears to show that there are some inconsistencies in the phenomenological 
arguments given above. To resolve these inconsistencies, we have attempted to examine the 
MHD equations using renormalization group analysis. 



For fluid turbulence Forster et al. || and Yakhot and Orszag || have applied dynamical 
RG procedure in which a forcing term with a power law distribution in wavenumber space 



is introduced. McComb ||, McComb and Shanmugasundaram ||, McComb and Watt [JIO 
and Zhou et al. |TI[ applied a self-consistent RG procedure that yields Kolmogorov's energy 
spectrum. For MHD turbulence, Fournier et al. |2IJ and Camargo and Tasso ||2T| have used 



RG procedure similar to that of Forster et al. || and Yakhot and Orszag ||. In all these 
schemes the averaging is done over the small scales (based on Wilson's approach in his 
Fourier space RG). Till date the RG methods applied to MHD turbulence do not find direct 
evidence of Kolmogorov-like power law in MHD turbulence. In a more recent work, Verma 
and Bhattacharjee [2^] have applied Kraichnan's DIA |||[] to MHD turbulence and obtained 



the Kolmogorov's constant for MHD, but in Verma and Bhattacharjee's work k~ 5 ^ 3 energy 
spectra was assumed, and an artificial cutoff was introduced for the self energy integral. 

In this paper we construct a self-consistent RG procedure similar to that used by Mc- 
Comb |J, McComb and Shanmugsundaram ||, McComb and Watt flIQ| , and Zhou et al. (11 



for fluid turbulence. However, one major difference is that we integrate the small wavenum- 
ber modes instead of large wavenumber mode integration used by earlier authors. In our 
procedure we obtain the effective mean magnetic field B (k) as we go from small wavenum- 
bers to large wavenumbers. At small wavenumbers the MHD equations are approximately 
linear. During the RG process, the effects of the nonlinear terms in the small wavenumber 
shells is translated to the modification of B (k) at larger wavenumbers. 

We postulate that the effective mean magnetic field is the magnetic field of the next- 
largest eddy contrary to the KD phenomenology where the effective mean magnetic at any 
scale is constant. To illustrate, for Alfven waves of wavenumber k, the effective magnetic 
field Bi(k) (after ith iteration of the RG procedure defined below) will be the magnetic field 
of the eddy of size k/10 or so. This argument is based on the physical intuition that for 
the scattering of the Alfven waves at a wavenumber k, the effects of the magnetic field of 
the next-largest eddy is much more than that of the external field. The mean magnetic 
field at the largest scale will simply convect the waves; the local inhomogeneities contribute 



to the scattering of waves which leads to turbulence (note that in WKB method, the local 
inhomogeneity of the medium determines the amplitude and the phase evolution). In our 
self consistent scheme we find that Bq appearing in the Kraichnan's or Dobrowolny et al.'s 
argument must be k dependent. The substitution of k dependent B (k) leads to k " 5//3 energy 
spectra, which is consistent with the solar wind observations and the simulation results. We 
will describe these ideas in more detail in the following section. 

The normalized cross helicity a c , defined as (E + — E~)/(E + + E~), and the Alfven 
ratio ta, defined as the ratio of fluid energy and magnetic energy, are important factors 
in MHD turbulence. For simplicity of the calculation, we have taken E + (k) = E~(k) and 
ta = 1. These conditions are met at many places in the solar wind and in other astrophysical 
plasmas. 



II. CALCULATION 



The MHD equation in the Fourier space is [T3[ 

{-iw =F i (B ■ k)) zf(k,u) = -iM ijm (k) J dpduj'zf(p, uj')z±(\l - P, <*> - u') (6) 

where 



-^ijm(k) — kjPi m {\i)] Pi m {\^) — 5i m ^ 2 , (7) 



Here we have ignored the viscous terms. The above equation will, in principle, yield 
an anisotropic energy spectra (different spectra along and perpendicular to B ). Solving 
anisotropic equations is quite complicated. Therefore, we modify the above equation to the 
following form to preserve isotropy: 

(-iu =F % (B k)) zf(k,uj) = -iM ijm (k) J dpdLu'zJ(p,Lj')z±(k - p,u - cu') (8) 

This equation can be thought of as an effective MHD equation in an isotropic random mean 
magnetic field. 
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In our RG procedure the wavenumber range (fco..fcjv) is divided logarithmically into N 
shells. The nth shell is (k n _i..k n ) where k n = s n k (s > 1). In the following discussion, firstly 
we carry out the elimination of the first shell (ko--kx) and obtain the modified MHD equation. 
We then proceed iteratively to eliminate higher shells and get a general expression for the 
modified MHD equation after elimination of nth shell. The details of the renormalization 
group operation is as follows: 



A. RG Procedure 

1. Decompose the modes into the modes to be eliminated (A; < ) and the modes to be 
retained {k > ). In the first iteration (&0..&1) = k K and {k\..kfi) = k > . Note that B (k) 
is the mean magnetic field before the elimination of the first shell. 

2. We rewrite the Eq. (H) for fc < and k > . The equation for ^ ±:> (k, t) modes is 



'—iu =F i (B k)) zf > (k, u) = -iM ijm (k) J dpdu' zj> (p^^z^ (k - p,u - J 



+ 



zf > (P> ^') z ti < (k - p, a; - lu') + zj< (p, J)z^ (k - p, u - <J) 
zJ K (p, w')^ (k - p, u - u)') 



+ 



(9) 



while the equation for z i < (k, t) modes can be obtained by interchanging < and > in 
the above equation. 

3. The terms given in the second and third brackets in the RHS of Eq. (|9D is calcu- 
lated perturbatively. We perform ensemble average over the first shell which is to be 
eliminated. We assume that zf K (k, t) has a gaussian distribution with zero mean. 
Hence, 



*T(M)) = o 



(10) 



and 



z a < (p, u/)C(q> tO) = ^ m ( P )^(p, «0*(P + q)*(w' + J' 



(11) 



where a, 6 = ±. Also, the triple order correlations ^^(k, oj)z^ l < (p, u')zf < (q, uo")j are 
zero. We keep only the nonvanishing terms to first order. For the relevant Feynmann 
diagrams, refer to Zhou et al. [II|. Taking r& — 1 and E + (k) = E~(k), the Eq. ([|) 
becomes 

{-iu^i{B Q k))zt > {k,u) = -iM ijm (k) J dpduo' [x^(p, oj')z±>(k - p,u - u/)] + 
H) 2 M ijm (k) L +q=k dqckj'MmMPj^G^p, u')C^<(q, uj - w')zt>(k, u) + 



{-if M ijm (k) J p+q=k dqdu'M mst (p)P js (q)G^(p, cu')C^<(q, u - u')z±>(k, uj 



(12) 



where G is the Green's function obtained from the equation 



G-\k,u) 



( 



\ 



-iu - ikB^ + (k) -ikB£~ (k) 
ikE>Q + (k) —iu + ikP>Q~ (k) , 



(13) 



In deriving Eq. ( O ) we have neglected the contribution of the triple nonlinearity 
zf > (k, u)z^ > {p, uj')zf > {q ) ui"). McComb, McComb and Shanmugsundaram, and Mc- 
Comb and Watt PHTDI have also ignored the triple nonlinearity for fluid turbulence. 



4. Since Ta — 1 and E + (k) = E (k), we find that Bq (k) = B + (k). We also assume 
that the correlation functions C ±± have the same frequency dependence as G ±=l= , i.e., 

C ±± (k) 



-iu ± ^ikB^ ± (k) 



(14) 



Note that C (k) = E ±± {k) / {Auk 2 ) in three dimensions. From dynamical scaling 
arguments 



c± = T kB^(k) 



±±t 



(15) 



After some manipulations the Eq. (^) becomes 
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-no =p i 



B (k) + 5B^(k)\ k) zt>{k,t) =f i5B± T (k)z? > {k,t) 
= M ijm (k) J dp [zf>(p, t)z±>(k - P , t) 



(16) 



where 



x 



«^(*) = -*Wdq(|$) 

a 2 (fc,p, Q ) (X Q ±± (p)+B** Qp)) -g 4 (fc J»,q)B+- (p) 



and 



2X (p)(fcB± ± (fc)+pX ±± (p)- O X± ± ( O )) 



^(fc) = -fc/ p4 ^_ k dq(|^)x 

«3 (fc J»,q)B+- (p)-ai (fc,p,q) (X p ±± (pHB** (p)) 



(17) 



2X ±± (p)(fcB± ± (fc)+pX± ± (p)-gX ±± (g)) 



(18) 



where 2k 2 ai {k,p,q) = Ai{k,p,q) andX ±± (fc) = ^(B^ik)) 2 - (5^ T2 (A;)) 2 . The terms 
Ai(k,p,q) are given in the Appendix of Leslie || as Bi(k,p,q). Since, -E + = E~ and 
r A = 1, it is clear that = SB Q -(k). Therefore, J B ++ (/c) = B Q -(k) = B (k) 

and X£ + (k) = X -(k) = X (k). 

Let us denote Bi(k) as the effective mean magnetic field after the elimination of the 
first shell. 



B 1 {k) = B Q {k) + 8B (k) 



(19) 



Similarly, 



B+-(k) = B+-(k) + 5B+-(k) 



(20) 



5. We keep eliminating the shells one after the other by the above procedure. After n + 1 
iterations we obtain 

B* 1 (k) = B*(k) + 6B»(k) (21) 

where the equations for 5B^ ± (k) and 5B^ T (k) are the same as the equations ( |TT| , |TBD 
except that the terms BQ b (k) and XQ b (k) are to be replaced by B^ b (k) andX^ b (A;) 
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respectively. Clearly B n+ i(k) is the effective mean magnetic field after the elimination 
of the (n + l)th shell. 

The set of RG equations to be solved are (pTT| , p~8| ) with B replaced by B n s, and (f2~T|). 



B. Solution of RG equations 



To solve the Eqs. (|l7l , |l"8D with B n s and (plf), we substitute the following forms for E(k) 
and B n {k) in the modified equations ([I7| . [l8| ) 



with k = k n+ ik' 
We use n+ = rr 
B* ab (k') that are 



(22) 



E{k) = KU 2 / 3 k- 5 / 3 
B ab {k n k') = K^U^k-^B^ik') 

k' > 1). We expect that B* n ab (k') is an universal function for large n. 
- IT due to symmetry. After the substitution we obtain the equations for 



^(fc') = -/ p , +qW dq'(3)x 

a2(k,p,q)(X„(sp)+B n (sp))-a4(k,p,q)B+~ (sp) 
2X n (sp')(k'B„(sk')+p'X n (sp')-q'X n (sq')) 



(23) 



^;+-(^) = -W =k ^q'(S^)x 

az(k,p,q)B+-(sp')-a 1 {k,p,q)(X n {sp l )+B n {sp')) ~\ 
2X n (sp')(k'B n (sk')+p'X n (sp')-q'X n (sq')) 

B* n %(k) = s 1/3 B* n ab (k) + s~ l 'HB™ b (k) (25) 

Now we need to solve these three equations self consistently. The integrals in the Eqs. 
([23|,[24]) is performed over a region 1/s < p',q' < 1 with the constraint that p'+q'= k'. We 
use Monte Carlo technique to solve the integral. Since the integrals are identically zero for 
k' > 2, the initial S*(^) = B l nitial for k\ < 2 and B*^) = B l nitial * (k'j2)- 1 / 3 for k! > 2. We 
take -Bo" - = 0. The Eqs. (|23||2^) are solved iteratively. We continue iterating the equations 
till B* n+l (k') w B* n (k'), that is, till the solution converges. For B l nitial = 1.0, the B f n s for 



9 



various n ranging from ... 3 is shown in Figure 1. Here the convergence is very fast, and 
after n = 3 — 4 iterations -£>*(&) converges to an universal function 

f(k') = 1.24 * A;'- 32 . 

From the above arguments, we have shown that B*(k') is approximately proportional to 
£./-i/3^ rj^g other parameter B^ + ~(k') remains close to zero. 

We infer from the above analysis that the mean magnetic field scales as k^ 1 ^ 3 , and 
the energy spectra scales as k~ 5 ^ 3 . Essentially, the scaling of B leads to k~ 5 ^ 3 energy 
(Kolmogorov-like) spectra in our scheme. We have calculated B*(k') for B l mtial = 1,2,10 
and found that for large n, B^(k') « 1.25-B£" itiaZ A; /-1 / 3 or 

B n {k) = 1.25B l nitial K 1/2 U 1/3 k~ 1/3 . (26) 



C. Calculation of K 

We can calculate the Kolmogorov's constant for MHD turbulence K by calculating the 
cascade rate II 0]. In MHD the cascade rates are 

n + (k) = IT (jfe) = - / dk'T(k') (27) 
Jo 

The numerical solution of the cascade rate integral yields [|] 

1 1/1 r>initial 

^^^ = 3-85 (28) 

From the above equation it is evident that the Kolmogorov's constant K is dependent on 
the mean magnetic field B % n%t%a \ in fact, K oc (£?* m * ia ') 2 / 3 . Clearly, an increase in the mean 
magnetic field leads to an increase in the Kolmogorov constant, which in turn will lead to 
a decrease in the cascade rate (cf. Eq. (|5|)). This result is consistent with the simulation 



results of Oughton [23]. However, a cautious remark is necessary here. We have considered 
the mean magnetic field to be isotropic; this isotropy assumption needs to be relaxed for 
studies of realistic situations. 
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III. CONCLUSIONS 



We obtain Kolmogorov-like energy spectrum in MHD turbulence in presence of arbitrary 
B by postulating that the effective B is scale dependent In our renormalization group 
scheme we find that the self consistent B n {k) is proportional to A; -1 / 3 and E{k) is proportional 
to /c~ 5 / 3 . This analysis has been worked out when E + = E~ and r A = 1. The generalization 
to arbitrary parameters is planned for future studies. 

In our methodology the averaging has been performed for small wavenumbers in contrast 
to earlier RG analysis of turbulence in which the higher wavenumbers were averaged out. 
Our scheme yields a power law solution for large wavenumber, and is independent of the 
small wavenumber forcing states. This is in agreement with the Kolmogorov's hypothesis 
which states that the energy spectrum of the intermediate scale is independent of the large- 
scale forcing. Any extension of our scheme to fluid turbulence in presence of large-scale shear 
etc. will yield interesting insights into the connection of energy spectrum with large-scale 
forcing. 

I thank V. Subrahmanyam, J. K. Bhattacharjee, and M. Barma for numerous useful 
discussions. 
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FIGURES 

Figure 1: B*(k') for n = ... 3. The line of best fit f(k') to B*{k') overlaps with B*. 
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